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Abstract 

Riemann-Hilbert techniques are used in the theory of completely integrable differ- 
ential equations to generate solutions that contain a free function which can be used 
at least in principle to solve initial or boundary value problems. The solution of a 
boundary value problem is thus reduced to the identification of the jump data of the 
Riemann-Hilbert problem from the boundary data. But even if this can be achieved, 
it is very difficult to get explicit solutions since the matrix Riemann-Hilbert problem 
is equivalent to an integral equation. In the case of the Ernst equation (the stationary 
axisymmetric Einstein equations in vacuum), it was shown in a previous work that the 
matrix problem is gauge equivalent to a scalar problem on a Riemann surface. If the 
jump data of the original problem are rational functions, this surface will be compact 
which makes it possible to give explicit solutions in terms of hyperelliptic theta func- 
tions. In the present work, we discuss Riemann-Hilbert problems on Riemann surfaces 
in the framework of fibre bundles. This makes it possible to treat the compact and the 
non-compact case in the same setting and to apply general existence theorems. 



PACS numbers: 04.20. Jb, 02.10.Rn, 02.30.Jr 



1 Introduction 



Riemann-Hilbert techniques provide a powerful tool if one wants to solve initial or boundary 
value problems for completely integrable differential equations. They are used to generate 
solutions with a prescribed singularity structure that contain a free function. In the case 
of a boundary value problem, this function has to be chosen in a way that the solution 
takes the prescribed values at the given boundary, and similarly for initial value problems. 
Whereas this does not pose any problems in principle, there is little hope in practice to get 
explicit solutions to boundary value problems in this way. The reason for this is that the 
free functions (the 'jump data' of the Riemann-Hilbert problem) enter the solutions of the 
differential equation as part of an integral equation from which the solutions have to be 
constructed. This implies that one has to solve this integral equation first on the boundary 
where the boundary data are prescribed in order to determine the jump data. In a second 
step, one has to solve the integral equation with the then fixed jump data in the region under 
consideration. In general, both steps cannot be done explicitly. 

In the case of the stationary axisymmetric Einstein equations, the situation is however 
different. The typical problem one has to consider there is the exterior of a relativistic star 
or a galaxy. Within these models, the matter leads to boundary conditions for the vacuum 
equations if a solution of the Einstein equations in the matter region is known. Special 
cases are two-dimensionally extended matter distributions where the field equations reduce 
to ordinary differential equations, e.g. disks that are discussed in astrophysics as models 
for certain galaxies. Thus one has to consider boundary value problems for the vacuum 
equations in the case of compact matter distributions where the boundary is the surface of 
the matter, and where the boundary data follow from the metric functions in the matter 
region. Since the vacuum equations are equivalent to a single complex differential equation, 
the Ernst equation [|IJ, which is completely integrable |2|, f|, f|], one can use Riemann-Hilbert 
techniques to solve the resulting boundary value problems. 

In a previous work || , it was shown that Riemann-Hilbert problems for the Ernst equa- 
tion with analytic jump functions are gauge equivalent to a scalar problem on a Riemann 
surface. In the case of rational jump data, this surface is compact which makes it possible to 
give explicit solutions to the Ernst equation in terms of hyperelliptic theta functions. Thus 
it is not necessary to consider integral equations in this case. The physical properties of the 
resulting class of solutions were discussed in || where it was shown that the solutions can 
have the expected regularity properties and asymptotic behaviour. In the present article, 
we discuss the Riemann-Hilbert problems on Riemann surfaces in the framework of fibre 
bundles which makes it possible to treat the case of compact and non compact Riemann 
surfaces within the same setting. Using a theorem of Rohrl 0, we obtain an existence proof 
for the solutions to the Riemann-Hilbert problems. In the case of non compact Riemann 
surfaces, the constructed bundles are trivial due to Grauert's theorem ||. 

The paper is organized as follows: In section 2, we recall that the Ernst equation can 
be treated as the integrability condition for an overdetermined linear differential system, for 
which we formulate the matrix Riemann-Hilbert problem. As an example we consider the 
scalar problem in the complex plane which can be solved with the help of the Plemelj formula 
0. The matrix problem for the Ernst equation is equivalent to an integral equation which 
cannot be solved explicitly in general. Explicit solutions can in general only be obtained if 
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the jump matrix is diagonal thus implying that the resulting solutions are static. In this 
particular case the Ernst equation reduces to the axisymmetric Laplace equation. By the 
help of gauge transformations of the linear system of the Ernst equation, we are able to 
transform the matrix problem to a scalar one on a Riemann surface in section 3. We discuss 
the relation between fibre bundles and Riemann-Hilbert problems on Riemann surfaces. 
These results are used to prove the existence for the solutions of the original problem. In the 
case of non compact Riemann surfaces, Grauert's theorem [BJ implies that the constructed 
bundles are trivial. On compact surfaces, we recover the explicit solutions in terms of theta 
functions. 



2 The Riemann-Hilbert Problem for the Ernst Equa- 
tion 

It is well known that the vacuum metric in the stationary and axisymmetric case can be 
written in the Weyl-Lewis-Papapetrou form (see [fLOp 



ds 2 = -e 2U (dt + ad0) 2 + e~ 2U (e 2k (dp 2 + d( 2 ) + p 2 d0 2 ) , (2.1) 

where p and ( are Weyl's canonical coordinates and d t and d<p are the two commuting (asymp- 
totically) timelike respectively spacelike Killing vectors. In this case, the field equations are 
equivalent to the Ernst equation for the potential / where / = e 2U + ib, and where the real 
function b is related to the metric functions via b z = — (i/ p)e 4U a^ z . Here the complex variable 
z stands for z = p + i(. With these settings, the Ernst equation reads 

1 2 

fzz 4" _|_ [fx fz) = j ~ jfzfz i (2-2) 

where a bar denotes complex conjugation in C. With a solution of the Ernst equation, the 
metric function U follows directly from the definition of the Ernst potential whereas a and 
k can be obtained from / via quadratures. 

The importance of the formulation of the field equations in terms of the Ernst potential 
arises from the fact that the Ernst equation is completely integrable, see 0, ^, |J. This means 
that the latter can be treated as the integrability condition of an overdetermined linear 
differential system that contains an additional complex parameter, the so called spectral 
parameter, that reflects an underlying symmetry of the Ernst equation. We use the linear 
system for the 2 x 2-matrix $ of Q, 

Here the spectral parameter K resides on a family of Riemann surfaces C(z, z) parametrized 
by the physical coordinates z and z and given by Pq{K) = (K — iz)(K + iz). A point on L 
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is denoted by P = (K, fx (K)) with K 6 C. The functions M and N depend only on z and 
z but not on K, and have the form 

M = lh' N = JTT (2 ' 5) 

where / is again the Ernst potential. 

To construct solutions to the Ernst equation with the help of the above linear system, 
one investigates the singularity structure of the matrices and with respect to 

the spectral parameter and infers a set of conditions for the matrix $ that satisfies the linear 
system ( |2.3j ) and (|2.4|) . This is done (see e.g. jTT)) in 



Theorem 2.1 Let $ be subject to the following conditions: 

I. 3>(P) is holomorphic and invertible at the branch points Pq = —\z and Pq such that the 
logarithmic derivative 3> 2 $ _1 diverges as (K + iz)? at Pq and as (K — 12)2 at Pq. 

II. All singularities of $ on C (poles, essential singularities, zeros of the determinant of 
$, branch cuts and branch points) are regular which means that the logarithmic derivatives 
^ z ^~ x and are holomorphic there. 

III. $ is subject to the reduction condition 

$(P CT ) = a^{P)a x (2.6) 

where a is the involution on £ that interchanges the sheets, and o~\ and are Pauli matrices. 

IV. The normalization and reality condition 

$(P = oo+)= ({ ]) (2.7) 



is fulfilled. Then the function f in $>7\) is a solution to the Ernst equation. 
This theorem has the following 

Corollary 2.2 Let $(P) be a matrix subject to the conditions of Theorem 2.1 and C(K) be 
a 2 x 2-matrix that only depends on K G C with the properties (the a, are scalar functions) 

C{K) = ai (K)l + a 2 (K)a l , 

ai(oo) = 1, a 2 (oo) = 0. (2.8) 

Then the matrix $'(P) = <&(P)C(K) also satisfies the conditions of Theorem 2.1 and 
$'(oo+) = $(oo+). 

In other words: matrices $ which are related through the multiplication from the right by 
a matrix C of the above form lead to the same Ernst potential though their singularity 
structure may be vastly different (the functions a, need not be holomorphic). Therefore this 
multiplication is called a gauge transformation. 

Theorem ^J] can be used to construct solutions to the Ernst equation by determining 
the structure and the singularities of $ in accordance with the conditions I to IV. In the 
present paper, we will concentrate on the Riemann-Hilbert problem for the Ernst equation 
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which can be formulated in the following form: Let T be a set of (orientable piecewise 
smooth) contours G C (k = 1, . . . , I) such that with P G T also P G V and P a G T. Let 
Gk{P) be matrices on with analytic components and nonvanishing determinant subject 
to the reality condition Qu(P) = Gu(P) for the diagonal elements, and Gij(P) = —Qij{P) 
for the offdiagonal elements. We define j(t, Fj) — 1 if t G Tj and otherwise and Q = 
Y}k=il{t,T k )Q k . Let G(P a ) = aiG(P)o~i. Both T and Q have to be independent of z, z. The 
matrix $ has to be everywhere regular except at the contour T where the boundary values 
on both sides of the contours (denoted by $±) are related via 

*_(P) = $ + (P)ft(P)| P6r4 . (2.9) 

It may be easily checked that a matrix $ constructed in this way satisfies the conditions 
of Theorem |2.1| . Furthermore, it can be seen from the Theorem that the only possible 
singularities of the Ernst potential can occur where the conditions are not satisfied, i.e. 
where $ cannot be normalized or where P$ coincides with one of the singularities of $, in 
our case the contour V. The latter makes the Riemann-Hilbert problem very useful if one 
wants to solve boundary value problems for the Ernst equation: choose the contour V in 
a way that Pq G T just corresponds to the contour in the meridian (z, z)-plane where the 
boundary values are prescribed. The Ernst potential will in general not be continuous at this 
contour, but its boundary values will be bounded. Notice however that the Ernst potential 
will not always be singular if P coincides with a singularity of $ since the latter may be 
e.g. pure gauge. Theorem [2.1| merely ensures that the solution will be regular at all other 
points. 

To solve Riemann-Hilbert problems, one basically uses the same methods as in the 
simplest case, the problem in the complex plane for a scalar function ip, see e.g. [T2\. If Tk 



is a simply connected closed smooth contour and G a nonzero Holder continuous function 
on in C, the function ip that is holomorphic except at the contour Tk where 

i). = ifj + G (2.10) 

is obviously given by the Cauchy integral, 



m) = F(K)e,^-J rii — j, (2.1i; 



where F(K) is an arbitrary holomorphic function, and where the principal value of the 
logarithm has to be taken. The well-known analytic properties of the Cauchy integral ensure 
that condition ( p.lOj ) is satisfied. Formula Q2.ll ) shows that the solution to a Riemann- 



Hilbert problem of the above form is only determined up to a holomorphic function. Since 
F is holomorphic, the solution will be uniquely determined (due to Liouville's theorem) by 
a normalization condition ijj(oo) = ipo for oo ^ T^- Uniqueness is lost if one allows for 
additional poles since F in ( |2.11| ) has then to be replaced by a meromorphic function. We 
note that the above conditions on the contour may be relaxed: T may consist of a set of 
piecewise smooth orientable contours which are not closed. The Plemelj formula, see e.g. 0, 
assures that formula fl2.11|) still gives the solution to (|2.10|) . A normalization condition will 



however only establish uniqueness of the solution if G — 1 at the endpoints of T, i.e. if the 
index of the problem (|2.10|) vanishes. 
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Riemann-Hilbert problems on the sphere C which occur in the case of the Ernst equation 
can be treated in much the same way as the problems in the complex plane. The basic build- 
ing block for the solutions is the differential of the third kind cLjk+k- {X) that corresponds 
to the differential dXj (X — K) in the complex plane, i.e. a differential that can be locally 
written as F(X, K)dX where F(X, K) is holomorphic except for X = K ± where the residue 
is ±1. If we make the ansatz 

$ " $0 = 27i L x( X )^k + k-(X), (2.12) 

where the 2 x 2-matrix % is given by % = SL=i t(^j ^k)Xk and where $ is holomorphic, we 
get for at the contour Y with the Plemelj formula $ ± = ±|x+ Jr x{X)duj K + K -(X). 
Thus the Riemann-Hilbert problem ([2.9D is equivalent to the integral equations (at V) 

\x(<3+l) + ^- I X (X)dw K+K -(X)(g -1) = . (2.13) 

For simplicity, we will only consider the case where the projection of the contour T into 
the complex plane has a simply connected component r^. In || we have shown that the 

general problem ( |2.9| ) is gauge equivalent to a problem with Q\ = ( °t ^ ] on T\ which is 



/I, 

the contour in the + sheet. This gauge transformation does not change the singularity 
structure of $ (i.e. $ will only be singular at V). The reduction condition III of Theorem 

2.1 implies that Q2 = ( g ^ J on the contour I^ in the — sheet. Because of the reality 



conditions for Q, this implies that solutions to the Ernst equation that follow from 
contain two real valued functions which correspond to a and (3. The reduction and reality 
properties of $ (see Theorem 2.1) make it possible to consider only one component of the 
matrix, e.g. $12 from which the Ernst potential follows as / = $i2(oo + ). With $12 = 
^° + /ri 2wfx)(x^K) -^(^' z -> z)dX we obtain the Ernst potential for given a and (3 where 
Z is the solution of the integral equation 



2 4vri Jn fJto(X)(X - K) v ' ' ' 4m J ri Ho(X)(X - K) 

(2-14) 

and where i/jq follows from the normalization condition $ 12 (oo ) = 1. 

Explicit solutions can in general only be obtained for diagonal Q, i.e. for (3 = 0. In this 
case we get with the above formulas f — f — e 2U with 

U = -—[ , ^ dK. (2.15) 

47T1 M l( K _Q2 + p 2 



Thus all solutions are real in this case which implies that they belong to the static Weyl 
class. Since the Ernst equation reduces to the axisymmetric Laplace equation for U if / is 
real, the function U in fl2.15| ) solves the Laplace equation. In fact one can show that the 
contour integral there is equivalent to the Poisson integral with a distributional density. It 
can also be directly seen from the expression ( |2.15| ) that the dependence on the physical 
coordinates p and ( enters through the branch points of the family of surfaces C 
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3 Riemann-Hilbert Problems on Riemann Surfaces 
and Vector Bundles 

In the previous section we recalled that matrix Riemann-Hilbert problems cannot be solved 
in general. Only for particular cases one can find an explicit form of the solution. In |H 
we have shown that in the context of the Ernst equation it is, however, possible to go one 
step further if one drops the condition that the gauge transformed matrix $' has the same 
singularity structure as the original matrix $ in (|2.9|) . It was furthermore shown there that 
the Riemann-Hilbert problem Q2.9Q is gauge equivalent to a problem with diagonal matrix 
Q' = diag(G, 1) on a two sheeted covering L of C, given by an equation of the form 

fi 2 (K) = F(K)H(K) , (3.1) 

where F(K) and H(K) are holomorphic functions. They follow from the jump matrix Q via 

F{K) (Qu - Q12 + G21 - Q22) (Q11 - Q12 - Q21 + Q22) 



H(K) (On + Qn — G21 — £22) (£11 + Q12 + Q21 + Q 



22 1 



(3.2) 



whereas the analytic jump function G can be expressed via the components of the original 
jump matrix Q by 



G + l 



(Cu-fc-fo + MOb + fc + fc + ft,) (3 . 3) 



G —I \ (Qn - Q12 + Q21 - Q22) (Qu + G12 - Q21 - Q22) 

By definition a Riemann surface is given by an equation of the form f(K,fi) = 0, where 
f(K, fx) is an analytic function of K and /x. If f(K, /1) is a polynomial in both variables one 
speaks of a compact Riemann surface. Therefore, for our purposes it is sufficient to solve a 
scalar Riemann-Hilbert problem on a four sheeted Riemann surface which is, dependent on 
the initial variables Qij, either compact or non compact (if the components of Q are rational 
functions, the surface will be compact). 

In the mentioned paper we restricted ourselves to the case of C being compact with genus 
g, where it is possible to give explicit solutions to the Riemann-Hilbert problem in terms of 
theta functions. Here we give a characterization of solutions to Riemann-Hilbert problems 
in terms of fibre bundle theory, which allows for a treatment of both the compact and non 
compact case. 

It is a well known fact, see |13|], that there is a relation between Riemann-Hilbert problems 
on Riemann surfaces and holomorphic vector bundles over them. To make the paper self 
consistent, we begin with a brief introduction into fibre bundle theory, see [H, [TJJ. 

Let us recall what a differentiable fibre bundle is. 

Definition 3.1 A differentiable fibre bundle (E, M, F,G,tt) is a 5-tuple consisting of 
I. A differentiable manifold E - the so called total space. 
i7. A differentiable manifold M - the so called base space. 
1/7. A differentiable manifold F - the fibre or typical fibre. 

IV. A surjective map 7r : E — > M , called projection. The inverse image vr _1 (p) = F p ~ F of 
p G M is called the fibre at p. 
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V. A Lie group G - the structure group - acting on F on the left. 

VI. An open covering {U} of M together with diffeomorphisms fa : Ui x F — > 7r _1 ([/j) ; such 
that n o fa{p, f) = p (p G Ui, f G F). We call fa the local trivialisation since (p^ 1 maps 
K~ l {Ui) onto the direct product Ui x F . 

VII. If we set fa, p (f) = fa(p, f) then fa tP : F — > F, p is a diffeomorphism. If Ui C\Uj 7^ we 
require that tij(p) = 4>^p4>j, P '■ F — > F be an element of the structure group G. Then (pi and 
4>j are related by a smooth map : Ui D Uj — > G as 

fa(pj) = fa(p,tij(p)f) ■ (3-4) 
We call the {t^} the transition functions. 

Let us say a few words about this definition. If we take a chart U of the base space M then 
ix^iJJi) is diffeomorphic to Ui x F with diffeomorphism 0" 1 : ix^iU^ — > Ui x F. If the 
intersection [/, fl C/j 7^ 0, there are two maps 0« and 0j on this intersection. Let u & E such 
that 7T (it) = p G C/j fl C/j. We then have 

^r 1 («) = (p,/o 

There is a map tjj : U^nUj —> G which relates and /j-: = Uj(p)fj. 

The transition functions can not be chosen arbitrarily in order to be the transition func- 
tions of a fibre bundle. They have to satisfy some consistency conditions: 

tuip) = identity map , (p G U) 

Uj(p) = tjiip)' 1 , (peUnUj) 
Uj(p)tjk(p) = Uk(p) , (p e U n n C/ fe ) . (3.5) 

A fibre bundle is called trivial, if all the transition functions can be taken to be identity 
maps. A trivial bundle is of the form E = M x F. 

The transition functions are so important because they contain all the information 
needed to construct a fibre bundle. Let us now show how starting from a 5-tuple 
(M, {U}, {Uj(p)}, F, G) a fibre bundle over M with typical fibre F can be constructed. 
Finding the bundle means finding unique E, n and fa from the above data. We define 

X = \J(Ui x F) , (3.6) 

i 

and introduce on X an equivalence relation ~ as follows. We say that (p, f) G Ui x F and 
(q, /') E Uj x F are equivalent, (p, /) ~ (q, f), if and only if p — q and /' = tij(p)f. The 
total space E of the fibre bundle is then defined 



E = X/~ . (3.7) 

We denote an element of E by the equivalence class [(p, /)]. The projection 7r : £ — > M is 
given by 

[(p,/)]-P , (3-8) 
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and the local trivialisation fa : U{ x F — > 7r is given by 

&:(p,/)-[(p,/)] , (3-9) 

with p E Ui and f E F. The above data reconstruct the bundle £7 uniquely. 

Let us now make the relation between Riemann-Hilbert problems on Riemann surfaces 
and vector bundles on them more explicit. 



3.1 The compact case 

First we will consider the case that the Riemann surface £ obtained as the double covering 
via the procedure described in |5J of the Riemann surface C is compact, i.e. we consider the 
Riemann-Hilbert problem 

^(P) = MP)G(P) (3.10) 
on a compact Riemann surface. Due to the fact that on such surfaces theta functions are 



the basic building blocks for the construction of meromorphic functions, see 16 1, we may 



express the solution to ( |3.10| ) in terms of these functions. But, in order to make contact with 



the case of L being non compact, we describe here how a solution to ( 3.10Q is connected with 



some line bundle over C. To this end we want to make use of the above result that a vector 
bundle over a manifold £ is completely determined by a triple (£, {Ui}, {Uj}). Let {Ui} be 
a covering of L and {(fti(P)} solutions to ( 3.10|) in the domains U, different from zero. In 



the domains U the original scalar Riemann-Hilbert problem is reduced to a problem on the 
complex plane C, which can be solved (see the discussion in section |2|). 

In other words, the functions <fii{P) are non vanishing on Ui and fulfil ( |3.10|) on the 



intersection r H Ui. If Y n U$ = then fa(P) is a holomorphic function in Ui. Let us now 
define some functions tfj(P) in U fl Uj by 



^=¥m< 1:111 



for P EUnUj. We have 



t-(p) - - M p M p ) _ ^(P) _ f+(p) , u) 

i.e. the functions tij(P) = tfj(P) = t^(P) do not jump at the contour T. It follows im- 
mediately from this definition that the tij fulfil the consistency conditions for transition 
functions (|3.5|) . Because the <pi{P) are non vanishing the functions tij(P) take values in C* 
(the complex numbers different from zero). Therefore, we may define a complex bundle with 
structure group C* and standard fibre C, i.e. a line bundle, Bq, over the compact Riemann 
surface t by the 5-tupel (£, {Ui}, {tij(P)}, C, C*). In other words: for £ being compact we 
may associate to a Riemann-Hilbert problem ( ft .ICQ a line bundle over this surface. 



3.2 The non compact case 

Let us now turn to the case t being non compact. Contrary to the compact case we do not 
have the calculus of theta functions associated to a Riemann surface at our disposal in order 
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to construct meromorphic functions. Nevertheless, we may perform a similar construction 
as in the compact case and construct a vector bundle for a Riemann-Hilbert problem given 
on this surface. The remarkable point is now that this vector bundle is, due to a theorem 
by Grauert |§, a trivial one. 

To be more precise, let £ be equipped with a covering Af = {Ui, i G /}, where / denotes 
some set of indices. Let us suppose that there exists a number N, the covering constant, 
such that any point P G £ belongs to no more than N domains of the covering, see [13]. We 



assume that the contour T is compact and closed, dividing £ into, in general, non compact 
domains. To simplify the discussion we are looking for solutions (j){P) with finite Dirichlet 
integral, i.e. 

J # Ad0 < oo , (3.13) 

C\A(T) 

where A(T) denotes some neighbourhood of T. Let <pi{P) be, as above, a non vanishing 
function on Ui and solving there the Riemann-Hilbert problem ( 3.1 0|) on the intersection 



r n Ui. If T fl C/j = then 4>i(P) is a holomorphic function in Ui. Similarly to the compact 
case we now define functions tfJP) in U fl Uj by 



Again, we have 



f - (P) _ Mg) _ W) _ h+(P) _ t+m f o Hl 



for P G UidUj, i.e. the functions Uj{P) = tfAP) = tZ(P) do not jump at T, analogously to 
the compact case. These functions also obey the consistency conditions ( j3.5|) and, therefore, 
we may associate to the Riemann-Hilbert problem (|3.10|) a vector bundle in the same manner 



as in the compact case. But, whereas in the compact case one does not know in advance 
what the global structure of the bundle space B G looks like, we have in the present case the 
following 

Theorem 3.2 (Grauert) Any complex line bundle over a non compact Riemann surface 
is trivial. 



Proof : The proof can be found in H. □ 

From this theorem it follows that for non compact £ the line bundle Bq associated the 
Riemann-Hilbert problem ( |3.10|) has the form 



5G 



£xC . (3.16) 



To conclude, we have shown that for £ being non compact (where one does not have an 
explicit solution of the Riemann-Hilbert problem in terms of theta functions) there is a 
simple geometric characterization of it in terms of fibre bundles over £. Due to its local 
properties the bundle approach allows to reduce the scalar Riemann-Hilbert problems on 
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C to problems on the complex plane C which can be explicitly solved. For non compact 
Riemann surfaces C the resulting total space is given as a direct product whereas in the 
compact case the total space is, in general, twisted in a non trivial manner. On the other 
hand, in this last case, an explicit formulation of the solution in terms of theta functions is 
possible. 
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